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Abstract 

In this note we describe a new method of counting the number of unordered factorizations 
of a natural number by means of a generating function and a recurrence relation arising from 
it, which improves an earlier result in this direction. 

1 Introduction 

Consider the natural number 18. It has 4 distinct "factorizations", namely 

18 = 2.3.3 = 2.9 = 3.6 = 18. 

Similarly, there are 9 ways of factorizing 36: 36 = 2.2.3.3 = 2.2.9 = 2.3.6 = 3.3.4 = 2.18 = 3.12 = 
4.9 = 6.6 = 36. Our problem is to find this number for any natural number n. Since we are not 
distinguishing between 2 • 9 and 9-2, such a factorization is called unordered. A partition of a 
natural number n is a representation of n as the sum of any number of positive integral parts, 
where the order of the parts is irrelevant. The number of such partitions of n is known as the 
partition function and is denoted by p{n). Likewise the function p*{n) denotes the number of 
ways of expressing n as a product of positive integers greater than 1, the order of the factors in 
the product being irrelevant. For convenience, p*{l) is assumed to be 1. Clearly p*{n) is the 
number of unordered factorizations of n. In 1983, Hughes and Shalit [B] obtained a bound for 
p*{n), namely, p*{n) ^ 2n^^ which was then improved to p*{n) ^ n by Matties and Dodd [7J in 
1986. By this time Canfield, Erdos and Pomerance P] modified a result of Oppenheim regarding 
the maximal order of p*{n). They obtained another bound for p*{n) and described an algorithm for 
it. An average estimate for p*{n) was given by Oppenheim [8] which was also proved independently 
by Szekeres and Turan [Q]. Finally, in 1991, Harris and Subbarao [5j came with a generating 
function and a recursion formula for p*{n). One may consider [1] and [3] for some problems 
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associated with p* (n) . For a list of values and computer programming one may consider the website: 
http : / / www .research . att . com / cgi-bin / access . cgi /as/nj as/ sequences / 1 (sequence no. A001055). 

In this note, we describe a new method for counting p*{n) and obtain a generating function for 
it which is followed by a recurrence relation that generalizes the one given by Harris and Subbarao 
[S] as well as the one given by Euler for p{n). The final recursion formula improves the one given in 
[5] as it contains less number of terms. It is important to note that we wish to develop an algebraic 
approach to the problem which might be helpful for other similar situations in future. Also we 
note some errors in describing an equivalent form of the recurrence relation in [5] . Throughout the 
note we denote set of all natural numbers, non-negative integers, integers and rational numbers by 
N, Z+, Z, Q respectively. 



2 Representation of numbers by polynomials 

Consider the monoid (N, •) of natural numbers under usual multiplication. For any natural number 
n, let S{n) be the submonoid of (N, •), generated by the set of prime factors of n, i.e., if the prime 
factorization of n is 

n = p'l'pT---PT, (2.1) 
where pi are distinct primes, G N for all i = 1, 2, . . . , /c, {k G N), then 

S{n) = {pYpI' ...pI" eN\n,eZ+, i = l,2,...,k}. 

We show that S{n) has an interesting algebraic structure. Define the partial ordering <' on S{n) 

by 

a <' b <J=^ a divides b. 

This ordering on S{n) is, in fact, a lattice ordering where a\J h = lcm(a, b) and a Ab = gcd(a, b) 
for all a,6 G S{n). Moreover this lattice is distributive and bounded below by 1. A monoid S is 
called a lattice- ordered semigroup if it has a lattice ordering that satisfies a{bV c) = abV ac and 
(6 V c)a = bay ca, for all o, 6, c G S. Now for all a,b,c G S{n), o{lcm(6, c)} = lcm(a6, ac). So we 
have the following proposition: 

Proposition 2.1. For any natural number n, {S{n), ■,<') is a lattice- ordered semigroup. 



Definition 2.2. Now corresponding to each natural number n we associate a polynomial in the 
polynomial semiring Z^[x] as 

f{x; n) = ni + n2X + n^x^ + • • • + n^x^^^ , 

where n has the prime factorization ()2.1|) . 
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Next we define a binary relation ^ on Z||[x] as follows: 
Let /(x) = oo + aix + a2X^ + • • • + amx"^ and g{x) = 69 + hix + 52X^ + • • • + hnx"^. Then 

/(x) ^ g{x) <^=^ m ^ n and ^ 6j for all i = 0, 1, 2, . . . , m. 

Clearly ^ is a partial ordering on Finally, let us denote the set of all polynomials in Z(j'[x] 

of degree less than k by Pk [x] . 

Theorem 2.3. (^^[2;],+,^) is a lattice- ordered semigroup which is isomorphic to {S{n), ■,<'), 
where n has the prime factorization 112. 

Proof. Let f{x) = ai + a2X + 032;^ + • • • + akX^^^ , g{x) = hi + b2X + b^x^ + • • • + b^x'^^^ £ 
Pk[x]. Then it is routine to verify that f V g = ci + C2X + csx^ + • • • + CkX^^"^ and f /\ g = 
di + d2X + d^x"^ + • • • + dkX^~^, where Cj = max{aj, 6,} and dj = min{aj, bi}. Thus ^) 
is a lattice. Obviously, is an abelian monoid where the identity element is the zero 

polynomial. Now choose f{x) = ^ ajX*~^, g{x) = ^ biX^~^ , h{x) = ^ CjX*~^ € Pk[x\. Then 

i=l i=l i=l 

k k 

f + {gV h) = {ai + max{6i,Ci})x*~^ = Yl (max{ai + bi,ai + Ci})x*-^ = {f + g) ^ {f + h). 

i=l i=l 

Therefore {Pk[x\,+, ^) is a lattice-ordered semigroup. 

Now define a map ifj : S{n) — y Pk[x\ by ipim) = /(x; m) for all m E S{n). That tp is bijective 

k k 

follows from Definition 12.21 Let mi = Y\ p\^\ m2 = Yl vt"^ ^ S{n). Then 

mi < 7712 

<^=^ mi divides m2 

<^=^ rii ^ r2i for each i = 1, 2, . . . , fc 

*^=^ rii + ri2X + ri3X^ H h rikx'^'^ ^ r2i + r22X + r23X^ H h r2kX^~^ 

f{x;mi) S f{x;m2) 
<^=^ '(/'(mi) ^ 'ip{m2). 

k k k k 

Also 77117712 = n Pi^'^''^*- Then /(x; 77717712) = X] (n« + ?'2i)^'~^ = Y '>'iiX^~^ + ^ r2jX*~"^ = 

i=l 1=1 i=l i=l 

f{x;mi) + f{x;m2). Thus ^{mim2) = ilj{mi) + ip{m2). Therefore 7/; is an isomorphism, i.e., 
(5'(n),-,<-) = (Pa;[x],+,^), as required. □ 

Definition 2.4. For any /(x) G ^([[2;], let p(/) denote the number of partitions of the polynomial 
/(x) in terms of addition of polynomials (not all distinct) in Z(|[x], where the order of addition is 
irrelevant. We assume p{0) = 1. 



For example, the distinct partitions of the polynomial 2 + x in Zo^[x] are 

2 + x = (l) + (l) + (x) 

= (1) + (1 + ^) 
= (2) + (x) 
= (2 + x) 

Note that /(x; 12) = 2+x and compare the above partitions with usual factorizations: 12 = 2.2.3 = 
2.(2.3) = 2^.3 = 12. 

Theorem 2.5. For any natural number n, p*{n) = p{f{x;n)). 

Proof. Let n € N and F{n) denote the set of all factors of n. Then {F{n), <') is a sublattice of 
S{n). On the other hand, definqj the set 

Sifix)) = {g{x)GZ+[x] \g{x)Sf{x)}. 

Then S{f{x)) is a sublattice of Pa,.[x], where /(x) = /(x; n) and n has the prime factorization (j2.1|) . 
By Theorem 12.31 it follows that the restriction of the map ip on F(n) is a lattice isomorphism from 
{F{n), <■) onto [S{f{x; n)), ^ ) . Indeed, let m = "PyJ)^ ■ ■ -Pk ^ F{n) for some ri, r2, . . . , G Zq. 
Since m is a factor of n, we have rj ^ nj for each i = 1,2, ... ,k. Hence 'ip{m) = f{x;m) = 
ri + r2X + r^x"^ + • • • + r^x*^"^ ^ ni + n2X + n^x"^ + • • • + nkX^~^ = /(x; n) which implies il>{m) G 
S{f{x;n)). Conversely, let g{x) G S{f{x;n)). Then deg5((x) ^ deg/(x) = k — 1. Let ^((x) = 
bi + 62X + b^x"^ + • • • + bkX^^^ for some 61, 621 • • • > ^fc G ZjJ'. Then 6j ^ for each z = 1, 2, . . . , A; 
and g{x) = g{x;p^iP^2---'Pk) — '^{P^iP^2---Pk)'' which implies ilj{F{n)) = S{f{x;n)). Thus 
we have {F{n),<') = (S'(/(x; n)), ^ ). Also since /(x;mim2) = /(x;mi) +/(x;m2) for all 
nil, 1712 £ there exists a one-to-one correspondence between factorizations n = mim2 . ■ ■ 

of n with partitions /(x; mi) + /(x; 7712) + • • • + /(x; mr) = /(x; mim2 . ■ ■ nrir) = f{x; n) of /(x; n). 
Thus we have p*{n) = p{f{x; n)). □ 

Corollary 2.6. Let p be a prime number and n G N. Then p*{p^) = p{n). 

Remark 2.7. It is clear that the value of p*{n) is independent of the particular primes involved 
in the prime factorization expression of n, i.e., if n has the prime factorization ()2.ip and m = 
12^ ^ where qi are distinct primes, then p*{m) = p*{n). Thus the polynomial f[x;n) 

in Definition 12.21 may be different for different arrangement of primes in the prime factorization 



^ S{f{x)) is called the section of f{x) in Zj[a::] 



expression of n. But the value of n)) remains same for each such arrangements. In particular, 

p{2 + x) = p*(2^.3) = 4 = p*{2.3^) = p(l + 2x). Indeed, the distinct partitions of the polynomial 
1 + 2x in [x] are 

l + 2x = (l) + (x) + (x) 
= {l + x) + {x) 
= (l) + (2x) 
= (l + 2x) 

More generally, p* {p^q) = 4 for any pair of distinct primes p, q. 

3 Generating function and recurrence relations 

Let n be a natural number. We know that the number of partitions, p{n) of n is given by the 
following classical generating function found by Euler: 

1 

(1-X)(1-;?)(1-X3) 



n-) = (3.1) 



oo ^ 



(3.2) 



n=l 

oo 



= n(i+E^"'j (3-3) 

1=1 n=l 

oo 

= l + ^p{n)x''. (3.4) 

n=l 

In the above equalities, since p.3|) provides all possible positive integral powers of x less than 
n (with all possible multiplicities), product of these terms produce the term as many times as n 
can be expressed as a sum of positive integers which is exactly the number of partitions of n, i.e., 
the term occurs p{n) times. So we get the coefficient p{n) of x" in (j3.4p . Similarly, if we wish to 
find the number of partitions of the polynomial /(x) in the polynomial semiring Zq'[x], we have to 
consider the product of summations which provides all possible polynomials in Zo^[x] less than /(x) 
(with all possible multiplicities) as indices. Thus we have the following formal generating function 
for p{f{x)): 

^(-) = n 

9 6 Z+[x]* 

oo 

= n ^^-^^ 

= 1+ Yl PU)ef^''\ (3.7) 
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where Z+[a;]* = Z+[x] \ {0}. 

Remark 3.1. The expressions (|3.5p - (|3.7|) are merely formal in the sense that for any particular 
f{x) € Z(|[a;], the coefficients of e^^^^ in either side are same. So the convergence problem does 
not arise here. However, if one insists on it, one may replace e by ei = - in which case ()3.5p -(f37 



are absolutely and uniformly convergent for all positive integral values of x. For example, consider 

•^(1) = n — huj- The product Y[ — is convergent if Yl e^^^'' is convergent. 

g& z+[x]* ^'"'^ gez+lx]* g€Z+[x]* 

(1) °° 

Now Yl ^1 X] which is absolutely and uniformly convergent for ei = ^ [1]. 

Now (13. 5p can be written in the form: 

= n — ^ — (3-8) 



g is primitive n = \ 

which is again by (|3.4p . 



oo 

= n (i+E^'H^"'^"^ ^^-^^ 

g is primitive 

So we have the following generating function for p{f{x)): 

oo 



g is primitive 

Using this we describe a method of calculating p{f{x)): 



oo 

p{f{x))= coefficient of e^(^) in JJ (l + ^p(n) e"^?^^)) . (3.11) 



g is primitive 

:)2 



Example 3.2. Let n = 12. Then n = 2 .3 and so the associated polynomial /(x; 12) = 2 + x. Now 
primitive polynomials less than or equal to 2 + x in are 1, x, 1 + x and 2 + x. So we have 

p(/(x;12)) = p(2 + x) 

= coefficient of e^"*"^ in 

(l+p(l)e+p(2)e2) (l+p(l)e-) (l+p(l)ei+") {l + p{l)e'+^) 

= coefficient of e^+'' in (1 + e + 2e2)(l + e^')(l + ei+^')(l + 6^+^) 

= coefficient of 6^+^ in (1 + e + 2e2)(l + + 6^+^ + 6^+^) 

= coefficient of 6^+^ in 1 + e + 26^ + + 2e^+^ + 46^+^ 

= 4. 
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Thus we get that p*{l2) = p{f{x; 12)) = 4. 

Remark 3.3. (i) Note that in each step of the above calculation, we omit the terms e'^^^^ whenever 
h{x) > f{x; 12), as these terms have no further contribution in forming the term e^^^'^'^\ 

(ii) By the process of calculating p{f{x)), we are getting all the values of p{g{x)) for all 
g{x) ^ f{x) in 'Lq[x\. For example, the above calculation gives us 

p{l) = 1, p{2) = 2, p{x) = 1, p{l + x) = 2. 
Next we wish to obtain a recurrence relation for p{f{x)). From (jS.Sp and (j3.7p . we get that 

(i+ E • n = ^- ^^-^^^ 

/ez+[x]* 9ez+[x]* 
Now taking formal derivative^ on both sides of (j3.12p we get, 

( E p^f) ^^^'^ • /'(^)) • n (1 - «'^'^) + 

(i+ E E {(-«^^^v(x)) n (l-e^^^^))}) = 

which implies 

oo 

(3.13) 

Then equating the coefficient of e-^^^^ of both sides we get 

k 

pU)f'{x) = ( E + E (EK/-^5)). (3.14) 

I Oo<f r=l 

where c(/) is the contenllf] of the polynomial f{x) and kg = max{r G N | f{x) > rg{x)}. Consid- 
ering p{0) = 1 and replacing rg by g we finally have 

p{f)f\x)= E Kg) 9\x) p{f - g)- (3.15) 

0<9g/ 

where \{g) = Yl, \ ■ Now since polynomials on both sides of ()3.15p are identical, we may equate 
coefficients of each power of x. So if f{x) = ni + n2X + nsx^ + • • • + n^x^"^ and 62(5) denotes the 



^The formal derivative of a polynomial h{x) = ao + aia; + a-zx^ + ■ ■ ■ + a^x'' G [a^] is defined by h' {x) — 

oi + 2a2X + 3a3X^ + ■ ■ ■ + kakx''~^. One can easily extend this definition for formal power series. The operator 

derivative is additive and follows Leibnitz' rule. It is routine to verify that the derivative of e*^^^' = e^'^^^h' {x). 
^i.e., gcd of all coefficients of the polynomial f{x). 
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coefficient of x in g{x) for each g € ^q[x] such that < g{x) ^ f{x), then equating constant terms 
in (|3.15p we have 

Pif) n2 = ^^9) b2{g) pU - g). (3.16) 

Also by Remark 12.71 one may rearrange coefficients of f{x) which does not change the value of 
Thus for each i = 1, 2, . . . , A;, we may rearrange coefficients of f{x) in such a way that will 
be the coefficient of x. Then we have 

p{f)n^ = Yl Kg)Hg)p{f-g), (3.17) 

where bi{g) denotes the coefficient of x*"-*^ in g{x) for each g € such that < g{x) ^ f{x). 

Therefore suitably multiplying the above relations by powers of x and adding we get 

Pif) fix) = Y ^(9) gix) Pif - g), (3.18) 

0<g^f 
ffez+ [x] 

Remark 3.4. (i) We first note that (|3.18p is a nice generalization of Euler's recurrence relation for 
pin), (n € N) which is given by 

n 

np(n) = ^(T(j)p(n - j), (3.19) 

i=i 

where o"(j) denotes the sum of all divisors of j and p(0) is assumed to be 1. Now (j3.19p is immedi- 
ately obtained from (j3.18p by putting /(x) = n (the constant polynomial), as we already have, by 
Corollarv 12.61 pin) = p*ip"'). Note that A(j)j = aij) for all j = 1, 2, . . . , n. 

(ii) The recurrence relation (|3.18p also generalizes the one obtained by Harris and Subbarao 
[5]. In Remark 2 of [5] (pg.477), they described an equivalent form of their recurrence relation as 
follows: 

Consider the vector d{n) = (ai, a2, ■ ■ ■ ctk) for the natural number n = qi^q^^ ■ ■ ■ %^ ■ Then 

p*(a) ||a|| = v\&~mm\V (3.20) 

0</3<Q 

where p*((3) = P*in), \\d\\ = Y[j=i'^jj K'^) — -| f i< x^.} -'^Z* ^^'^ (3 < a means that f3j < aj for 

We first note that the limit under sum in (j3.20p should be < /3 < a as A(0) is not defined. 
Secondly, ||a|| will be Ylj=i^j^ verify p*(18) = 5 by (j3.20p . which is wrong. Finally, 



8 



while defining the ordering for vectors one has to use < instead of <. However keeping aside these 
printing mistakes, the correct version of (|3.2U|) is given by 

P*{S)\\d\\= p*ia-P)\mm, (3.21) 

6<0<a 

k 

where ||a|| = ^ aj. Now (|3.21|) immediately follows from (|3.18|) by putting x = 1 (or, adding 
i=i 

equations (|3.17p for i = 1, 2, . . . , A;). 

Now for any n E N with prime factorization (|2.1|) . we define c(n) = gcd {ni, n2, . . . n^} and 
K''^) — Then equating constant terms in (j3.18p we get 

i|c(n) 

Pif) ni = J2 ^(9) bi{g) p{f - g), (3.22) 

9ez+[x] 

where bi{g) denotes the constant term of the polynomial g{x) for each g S ^o^[a;] such that < 
g{x) ^ f{x). This implies 

p*{n)m = nX{d)p*{^), (3.23) 

d|n, pi |d 

where the summation runs over d = PYP2' ■ ■ ■ P^k < ri ^ ni and ^ ^ rij, j = 2, 3, . . . , k. 
Thus we have 




(3.24) 



Remark 3.5. We note that (j3.24p improves the recurrence relation obtained by Harris and Subbarao 
[S] as it contains less terms. In fact, the number of terms in (|3.24p is (n2 + ^){n^ + 1) . . . (n^ + 1) 
less than that of (j3.2ip for n = p^-^p^^ ■ ■ ■ p'^'' ■ Thus, in view of Remark 12.71 it is advisable to 
arrange prime factors of n in such a way that ni should be minimum among all rij's for quicker 
computation. 

We summarize the above results in the following: 

Theorem 3.6. Let n = p^^p^^ ■ ■ ■p^'^ be the prime factorization of a natural number n, where pi 
are distinct primes, nj G N for all i = 1,2, ... ,k, {k (z N), with ui = minjnj | i = 1,2, . . . , k}. 
Then 

p{f) f{x;n) = Y ^^9) g{x) p{f-g), (3.25) 

OKg^f 
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where f{x;n) = ni + n2X + n^^x^ + • • • + UkX^ ^, X{g) = X] '^(S') ^ein^ the content of the 
polynomial g{x) for each g G ^o^[a^] such that < g{x) ^ f{x). In particular, 



p\d / 



where A(m) = X] l' c{m) = gcd {i,r2, ■ ■ ■ r^} form = P1P2 ■ ■ ■ P^k ^ < i ^ ni and Q ^ ^ 

j|c(m) 

raj, j = 2,3, . . . , A;. 

A simple formula is obtained the case of ni = 1. In that case, c(m) = 1 for all m\n with pi\m 
and So A(m) = 1 for all such m. Thus we have 

Corollary 3.7. Let n he a natural number and p he a prime number such that p /( n. Then 

p*{np)=Y,P*{d)- (3.26) 

In particular, for any two distinct prime numbers p, q and for any natural number n, 

n 

v\v<t) = Y.V^')^ (3-27) 

1=0 

Example 3.8. Let n = 72. Then n = 2^.3^ = 3^.2^. Now the divisors of 2^ are 1,2,4,8. So we 
have 

p*(72) • 2 



= {p*(24) +p*(12) +p*(6) +p*(3)} + 

2 • p*{8) + A(32 • 2) p*(4) + A(32 • 2^) p*(2) + X{3'^ ■ 2^) p*(l)}. 

Now p*{2) = p*{3) = 1, p*(4) = p*(22) = p{2) = 2, p*(8) = p*{2^) = p{3) = 3 by Corollary[Ml 
Again p*(6) = p*(3 • 2) = 1 +p(l) = 2, p*(12) =p*{3- 2^) = 2 + p{2) = 4, p*(24) = p*(3 • 23) = 
4 + p(3) = 7 by ([3211) ■ Thus we have 

p*(72) • 2 = {7 + 4 + 2 + 1} + 2 • {(1 + ^) • 3 + 1 ■ 2 + (1 + i) • 1 + 1 • 1} = 14 + 18 = 32. 

Hence p*{72) = 16. Indeed 16 factorizations of 72 are 

72 = 2.2.2.3.3 = 2.2.2.9 = 2.2.3.6 = 2.2.18 = 2.3.3.4 = 2.3.12 = 2.4.9 
= 2.6.6 = 2.36 = 3.3.8 = 3.4.6 = 3.24 = 4.18 = 6.12 = 8.9 = 72. 
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